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Abstract 

High energy multiparticle photoproduction off nuclear targets is studied. The photon is 
assumed to interact in direct scattering processes or as a resolved gg-state according to the 
Generalized Vector Dominance Model. In the description of resolved interactions multiple soft 
^ , and hard processes in each gg-nucleon interaction are taken into consideration. The model, 
j> [ formulated within the framework of the two-component Dual Parton Model, is shown to agree 
' w ith hadron production data from real and weakly virtual photon-nucleus interactions. Dif- 
5_i ■ ferences between multiparticle production in photon-nucleus and hadron-nucleus interactions 



and features of jet production in photon-nucleus collisions at HERA-energies are discussed. 
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1 Introduction 



Many features of high energy photoproduction off hadrons are well understood within the 
QCD-improved parton model: the photon may couple directly to a parton of the hadron [direct 
processes) or it may enter the scattering process as a hadronic fluctuation {resolved processes) 
(recent reviews are given in [[l], |2|). Since the lifetimes of these hadronic fluctuations are 
typically long enough to develop properties of ordinary hadrons they can interact in soft and 
hard resolved processes |3], ||. The classification into direct and resolved photon interactions 
has been confirmed by many experiments |5|, 0- Moreover, it has been shown that both 
classes of processes exhibit different properties concerning multiparticle production || £|, [|, [|, 
p. Here, we extend a previous study of photon-nucleus cross sections and their high energy 
shadowing behavior to particle production. This is done by applying the two-component 
Dual Parton Model (DPM) [TT|, [12], [TJ, [TJ, |15[ to photoproduction off nuclei. 

The two-component DPM is based on the DPM describing soft particle production (for a 
HI) and treats high-p^ processes using lowest order perturbative QCD. This model 



review see 



has been proven to be very reliable in describing the main features of hadron-hadron [|T|, 
|13|| , photon-hadron ]14|, |l5j, photon-photon [[Tj|, hadron-nucleus |17 |, and nucleus-nucleus 



collisions ||17|| . Here, the two-component DPM is applied to photon-nucleus collisions for the 
first time. The two channels which are of particular importance as a basis for our present 
study are hadron-nucleus and photon-hadron collisions. 

In hadron-nucleus collisions the multiple interaction process of the hadron with target nu- 
cleons which explains the "shadowing" -effect, is understood in the framework of the Gribov- 



Glauber approximation |1| [19|, p0[ . The two-component DPM incorporates the Gribov- 
Glauber approximation and treats each hadron-nucleon interaction by pomeron and reggeon 
exchanges. Dtunuc 2.0, a Monte Carlo (MC) realization of the two-component DPM for 



hadron-nucleus and nucleus-nucleus collisions, includes the Glauber-formalism the Pho 



jet model JTJ], [15| for each hadron-nucleon or nucleon-nucleon collision, a formation zone 



intranuclear cascade model, and models for spectator deexcitation and disintegration [22, E? 



As it will be shown further below, results on particle production in hadron-nucleus collisions 
obtained with Dtunuc 2.0 are in good agreement with data. 

The description of photon-hadron interactions within the two-component DPM takes the 
above mentioned dual nature of the photon into consideration. In analogy to hadron-hadron 
interactions, the resolved photon interactions are described in terms of multiple soft and hard 
pomeron and soft reggeon exchanges. Parton jet and minijet production in resolved processes 
as well as direct processes are calculated in lowest order perturbative QCD. The MC realization 
Phojet has been compared in several applications to photoproduction data from the HERA 
collider and a reasonable agreement was found || ^4], ^5|. 

Using all informations from hadron-nucleus and photon-hadron collisions it can be ex- 
pected that the two-component DPM also gives a reasonable description of multiparticle 
photoproduction off nuclei up to very high energies. Such a model might be useful for the 
study of general properties of photon-nucleus interactions as well as for the calculation of 
theoretical predictions (for example at HERA energies), cosmic ray cascade simulations, or 
shielding problems at TeV-lepton colliders. 

In Sect. 2 we summarize the basic ideas of the description of hadron-nucleon and photon- 
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nucleon interactions within the two- component DPM. Furthermore, its application to photon- 
nucleus interactions is presented. In Sect. 3 we demonstrate that the model is able to describe 
characteristic features of measured particle production in photon-nucleus interactions. The 
nuclear dependence of multiparticle photoproduction off nuclei and direct as well as resolved 
jet production are studied. Predictions for HERA-energies are given. Finally, in Sect. 4 we 
summarize our results. 



2 The two-component Dual Parton Model 

2.1 Photon— nucleon interactions 

In the following, we give a summary of the basic ideas of the two-component DPM with em- 
phasis to the description of photon-nucleon interactions as implemented in the MC generator 
Phojet. More detailed discussions can be found in Refs. [0, |I~5|j . 



The physical photon state is approximated as a superposition of a bare photon and of 
virtual hadronic states having the same quantum numbers as the photon. The bare photon 
may interact in direct processes. This direct contribution is estimated by lowest order pertur- 
bative QCD. The interactions of the hadronic fluctuations of the photon are called resolved 
processes and are described within the framework of the two-component DPM in terms of 
reggeon and pomeron exchanges. We distinguish soft and hard resolved processes according 
to the transverse momenta of the intermediate states of the pomeron exchange graph: By 
definition, scattering processes resulting in at least one parton having a transverse momen- 
tum larger than a cutoff momentum p™ toff = 3 GeV/c are called hard interactions. All other 
processes are classified as soft interactions. 

The amplitude describing the scattering of the hadronic fluctuation V with a nucleon iV 



is parametrized using a two-channel eikonal formalism [[LJ, pT3[] . In impact parameter repre- 
sentation, the amplitude reads [14 



av N (s, b) = - (l - e-*( s ' 6 >) , (1) 

with the eikonal function 

X(s, b) = xs(s, b) + xh(s, b) + xd(s, b) + xc(s, b). (2) 

Here, the Xi s denote the contributions from the different Born graph amplitudes: soft pome- 
ron and reggeon (S), hard pomeron (H), triple- and loop-pomeron (D), and double-pomeron 
amplitudes (C). The Born graph amplitudes of the soft processes are calculated using Gribov's 
reggeon field theory |fj6|, [L9| . Assuming universality of soft interactions, the free parameters of 
the model (pomeron couplings, pomeron intercept, triple-pomeron coupling, and slope param- 



eters) have been determined by a global fit to pp, pp, and jp cross section and slope data |T4 
However this universality cannot be applied to photon interactions involving large transverse 
momenta. Lowest order perturbative QCD is used to calculate the resolved and direct hard 
photon cross sections 0, j27j] . 
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The optical theorem relates the total cross section (i.e. elastic scattering, diffraction dis- 
sociation, and nondiffr active inelastic scattering) to the discontinuity of the amplitude ([I]) 
taken at vanishing momentum transfer. At high energies, the discontinuity can be expressed 
as a sum over terms which correspond to graphs with a certain number of cut pomerons 
(Abramovski-Gribov-Kancheli cutting rules [EH]). Therefore, one obtains for the exclusive 
cross section for k c cut soft pomerons, l c cut hard pomerons, m c cut triple- and loop-pomeron 
graphs, and n c cut double-pomeron graphs 

a{kalMS , b) = (^(^(^« e -«,>, (3 ) 

with 

J d 2 b a( y k c ,l c ,m c ,n c ,s,b) = a\; t N -a^ N . (4) 

k c +l c +m c +n c =l 

Both resolved soft and hard interactions are unitarized, i.e. a single gg-nucleon scattering may 
be built up of several soft and hard interactions. In contrast to resolved processes, absorptive 
corrections to direct photon interactions are suppressed by one order of the fine structure 
constant a em and are therefore neglected. 

In order to relate the cross sections given in Eq. @ to final state parton configurations, we 
use the equivalence between the reggeon and pomeron exchange amplitudes and certain color 
flow topologies p9 |. In Fig. |l] the color flows of single reggeon (a) and single pomeron exchange 
processes (c) are shown together with the corresponding cuts. Whereas a cut reggeon yields 
one color-field chain (string) (Fig. ]l]b), a cut pomeron results in two strings (Fig. |l]d) which 
are assumed to fragment independently into hadrons. In the large N c limit of QCD (with N c 
being the number of colors) the same color flow picture of the pomeron is also found in hard 
interactions. 

For pomeron cuts involving a hard scattering, the parton kinematics as well as the flavors 
and colors are sampled according to the QCD-improved Parton Model using leading-order 
matrix elements. We use the GRV parametrizations of the parton distribution functions 
(PDF) of the photon [[30| and the nucleon Both, initial and final state parton showers are 
treated. For the latter the Monte Carlo (MC) realization as implemented in Jetset |3j| ^ 
is applied. 

In case of pomeron cuts without large momentum transfer, we assume the partonic inter- 
pretation of the DPM: hadronic fluctuations of the photon are split into a quark-antiquark 
pair whereas baryons are approximated by a quark-diquark pair. The longitudinal momentum 
fractions of the partons are given by Regge asymptotics [[34], [35], [36], [37]. The valence quark 
momentum fraction x inside a nucleon is sampled according to p(x) ~ (1 — x) 15 /y / x whereas 
the momentum fraction of a valence quark inside a hadronic fluctuation of a photon is ob- 
tained from p(x) ~ 1/ Jx(l — x). For multiple interaction events, the sea quark momenta are 
sampled using a p(x) ~ 1/x behavior. The transverse momenta of the soft partons are not 
predicted by the DPM. Here, we assume an exponential distribution d 2 N s /dp 2 ± ~ exp(— f3p±). 
The energy-dependent slope parameter f3 is obtained requiring a smooth transition between 
the transverse momentum distributions of the soft constituents and the hard scattered partons. 

For single diffractive and central diffractive processes, the parton configurations are gener- 
ated using the ideas described above applied to pomeron-photon/hadron/pomeron scattering 
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processes ||38|| . Hence, a diffractive triple-pomeron or loop-pomeron cut can also involve hard 
scattering subprocesses resulting in a rapidity gap event with jets. 



All strings are hadronized using Jetset 7.3 [32, 33 . 

Finally it should be mentioned that the model is limited to photon virtualities Q 2 
with Q 2 < 9 GeV 2 [0. 



< s 



2.2 The event generator DTUNUC 2.0 for photon— nucleus colli- 



sions 



The MC event generator Dtunuc 2.0 for hadron-, photon-, and nucleus-nucleus collisions 
is based on its previous version p9l EUI incorporating the following new features: 



Dtunuc is extended to the description of photoproduction off nuclei. 

Particle production in each photon-, hadron-, and nucleon-nucleon interaction is based 
on the MC realization Phojet H, [15[ of the two-component DPM which includes 



multiple soft and hard pomeron exchanges and diffraction. 

A formation zone intranuclear cascade is treated in both, the projectile and target specta- 
tors including evaporation of nucleons and light nuclei from excited spectators, spectator 



deexcitation via photon emission, and fragmentation of light spectator nuclei [22, 23 



The generation of a photon-nucleusQ particle production event proceeds as follows: The 
treatment of any interaction starts with sampling of an initial spatial configuration accord- 
ing to nuclear densities and of the Fermi-momenta of the nucleons. The number of nucleons 
involved in the scattering process with the photon is obtained according to the Glauber ap- 
proximation. We use the MC algorithm of ||21|| , here extended to photon projectiles. The 
inelastic cross section for the scattering of a photon with virtuality Q 2 with a nucleus A at a 
squared photon-nucleon cm. energy s and at impact parameter b reads |10 



r M i ( 
a lA (s } Q 2 , b) = 4vr« em / dM 2 D(M 2 ) 



a VA (s,Q 2 ,M 2 ,b) = J Hd 3 rj p A (fA ( 1 



J^) (l + «£W.,«W)(5) 



i=l 



(6) 



Corresponding to the assumptions of the Generalized Vector Dominance Model (GVDM) (see 
for example [[3], |4]J and references therein) we integrate in (|5|) over the masses M of the 
hadronic fluctuations of the photon. D(M 2 ) denotes the density of hadronic states per unit 
mass-squared interval, e is the ratio between the fluxes of longitudinally and transversally 
polarized photons. Eq.((0) relates the inelastic cross section ay a for the interaction of the 



4 Here, we restrict our discussion to photon projectiles. The treatment of hadron-nucleus collisions is 
similar; the discussion of nucleus nucleus collisions is not within the scope of this paper 
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hadronic fluctuation with the nucleus to the amplitude V for the interaction of the hadronic 
fluctuation with a nucleon 



a yw (s,Q 2 ,M 2 



47nB(s, Q 2 , M 2 ) r/ r \2B(s,Q 2 ,M 2 ) 

&vn(s,Q 2 ) 
M 2 + Q 2 + C 2 



We refer to |T(J for the calculation of the M 2 -independent quantity 5"yiv, the slope 5, and 
the parameters entering the expressions. The integration over the coordinate space in (|6]) is 
performed using one-particle Woods-Saxon density distributions pa [fH]- With the impact 
parameter b and the mass M of the hadronic fluctuation which are sampled according to 
Eqs.([5],|6D, we obtain for the probabilities pi that an inelastic interaction between the projectile 
and nucleon i takes place |I| 

Pi^rfy + r*(b i ) - rfyrik) . (9) 

Note that bi is the impact parameter for the interaction of the hadronic fluctuation with the 
nucleon % for a fixed spatial configuration of nucleons. 

As it has been discussed in Sect. 3 of Ref. fllU| , hadronic fluctuations with M 2 > (2p™ toff ) 2 
predominantly interact in point-like interactions^, i.e. in direct interactions and in interactions 
which are characterized by the anomalous component of the photon PDF. In both cases, the 
virtuality of the qq-system allows to calculate the photon-gg coupling perturbatively. We 
denote the corresponding cross sections with cr^y and aOffi and refer to [JITJ for their calculation 
within lowest order perturbative QCD. Here we shall again consider the extreme case that in 
these processes the hadronic fluctuation of the photon interacts with only one target nucleon. 
Therefore, in the fraction A-(a^ + a^)/a^X of ai l events only one hard resolved (representing 
the anomalous component) or one direct photon interaction is sampled. 

Furthermore, as for the calculation of cross sections p0[ , the coherence length of the 
photon is taken into consideration which effectively leads to a suppression of the Glauber- 
cascade at low energies. Particle production in each inelastic interaction of the photon or 
its hadronic fluctuation with a nucleon is treated by the MC realization Phojet of the two- 
component DPM (see Sect . [Op . The photon-nucleon interactions are followed by a formation 
zone intranuclear cascade in the target spectator and by subsequent evaporation processes of 
nucleons and light nuclei as well as by spectator deexcitation, and by the fragmentation of 



light spectator nuclei [22, 23 



3 Particle production 

For studying photoproduction off nuclei we assume the logical sequence to be as follows: 
As mentioned earlier soft hadronic interactions exhibit universal features irrespective of the 

5 Note that in our definition of "point- like photon interactions" [|l0| interactions of low-mass qq- fluctuations 
with only one target nucleon are not included 
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nature of the colliding particles. This fact was already emphasized by Engel and Ranft in ] |T5|j 
where the description of high energy hadron-hadron and photon-hadron interactions based 
on the two-component DPM was extended to photon-photon interactions. Here, we want to 
proceed in an analogous way. We start from the description of photon-hadron and hadron- 
nucleus interactions and show that the model may also provide a reasonable description of 
photon-nucleus interactions. 

Photon-hadron collisions have been discussed in the framework of the two-component 



DPM elsewhere [|14l, 1131. Furthermore, also numerous studies of hadron-nucleus collisions 



within this model exist (see for instance |39], f4(J, [17j and references therein). Of course for 



the latter interaction channel it has to be verified that the new version of Dtunuc is able to 
describe the data. Since this has not yet been demonstrated we give a few examples further 
below. 

Based on the reasonable description of these two interaction channels the model is then 
applied to multiparticle photoproduction off nuclei. This application may serve as a severe 
test of the model since it carries no further freedom. 



3.1 Hadron-nucleus interactions 

In Table || average shower particle multiplicities calculated for interactions of positively charged 
pions and protons with different target nuclei at 50, 100, and 150 GeV are given. Here, shower 
particles are defined as singly charged particles with Lorentz-/? values exceeding 0.7. These 
results are compared to data taken from Ref. f43fl . For all configurations there is a reasonable 



agreement between the calculated and measured values. 

The distribution of the multiplicity of shower particles produced in interactions of 525 GeV 
pions in nuclear emulsions is plotted in Fig. ^]a. Model results shown as histogram are com- 
pared to data f44j[ . As this plot demonstrates also the distribution of the shower particle 



multiplicity is reproduced by the model. 



Calculated and measured [44] pseudorapidity distributions of shower particles are shown 
in Fig. ^3 again for pion-emulsion interactions at 525 GeV and, additionally, for 60 GeV. 
The invariant n° cross section as a function of the transverse momentum in proton-gold 



collisions at 200 GeV has been measured by the WA80-Collab. p5| . In Fig. ^|we compare our 



calculations to these data and find a reasonable agreement. Corresponding to the kinematic 
cuts of the experiment the pseudorapidity range 1.5 < rj < 2.1 is considered only. 



3.2 Photon— nucleus interactions 
3.2.1 General properties of the model 

Before comparing model results to data from lepton-nucleus interactions let us first out- 
line differences between hadron- and photon-nucleus collisions with respect to properties of 
the Glauber-cascade and multiplicities at fixed projectile energy. In lepton-nucleus interac- 
tions the projectile photon has varying energies and virtualities and a direct comparison with 
hadron-nucleus collisions would be less conclusive. 
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In Fig. |] we present average numbers of target nucleons interacting with the projectile {v t ) 
and multiplicities of shower (N s ) and heavy particles (N h ) (charged particles with (3 < 0.7) 
calculated for pion-, real photon-, and weakly virtual photon-copper interactions at energies 
of the projectile in the nucleus rest system between 10 GeV and 100 TeV. Differences between 
the z/ r values for pion and photon projectiles at fixed energy arise from differences between 
the pion-nucleon cross section and the (averaged over all masses M) gg-nucleon cross section 
(Tvn (Eq-dD) which enter the calculations [10|. Since the latter is slightly smaller, the real 



photon interacts with less nucleons than the pion. 

Furthermore, (z/ 4 ) decreases with the photon virtuality at fixed energy. There are two main 
effects being responsible for this behavior: 

(i) Due to the l/(M 2 + Q 2 )-behavior of ay at this cross section decreases with Q 2 at fixed M 2 
and fixed energy. Moreover, at large photon virtualities interactions of resolved photons 
with large masses M 2 become more important leading to a further decrease of the cross 
section <tvn and therefore (u t ). 

(ii) The fraction of events with point-like photon interactions [y t = 1) increases with rising 
photon virtuality. 

As an example, this is shown in Fig. [5]a for photon-copper interactions at a photon-nucleon 
cm. energy of 150 GeV, an energy which could be available in the future for photon-nucleus 
collisions at HERA jlfj . In our model, the average heavy particle multiplicity is approximately 
proportional to u t irrespectively of the nature of the projectile. This dependence is shown in 
Fig. |5|b for protons, pions, real and weakly virtual photons (Q 2 = 2 GeV 2 ). It reflects the fact 
that the number of heavy particles depends only on the number of nucleons "knocked out" 
of the target by the projectile and by subsequent intranuclear cascade processes and not on 



particular properties of the projectile-nucleon interactions [22, 



3.2.2 Average multiplicities 

Average shower particle multiplicities in interactions of 150 GeV muons in emulsions are 
compared to data |^7| in Fig. ^. Corresponding to the measurements, from the calculated 



events only those with more than two heavy particles in the final state have been taken 
into consideration. The flux of virtual photons is sampled according to the EPA and to the 
Q 2 -dependence of the photon-nucleus cross sections [10]. The multiplicities are shown as 
function of the inverse of the Bjorken-x. In this representation the average photon virtuality 
is decreasing from about 8 GeV 2 in the lowest bin (1/xBj < 10) to about 0.9 GeV 2 in the 



highest bin | 47| . Our results agree with the data in the kinematic region of photoproduction 
(Q 2 ~ 4 GeV 2 , 1/xbj > 25) but underestimate the data for higher virtualities. This fact 
may indicate that the PHOJET realization of the two-component DPM for photon-nucleon 
interactions fails in describing particle production in interactions of photons with relatively 
high virtualities at low energies. However, whereas the calculations are based in total on 50000 
events one has to note that the data include about 17 events in each 1/xBj-bin |[47|| . Taking 
into account that these 17 events involve interactions with six different targets (emulsion 
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components) at different collision energies and photon virtualities one might conclude that 
the statistical significance of the data is limited. 

Finally, we compare our results on multiplicities in muon-deuterium and muon-xenon in- 
teractions to data of the E665-Collab. f48|. The experiments were performed with a 490 GeV 



positive muon beam. The kinematic region under investigation is Q 2 > 1 GeV 2 , 8 < W < 
30 GeV, XBj > 0.002, and 0.1 < v/E^ < 0.85 (with v being the photon energy in the labora- 
tory). Furthermore, only charged particles with momenta p > 200 MeV/c in the laboratory 
frame are considered. We note that although the distribution of photons in leptons decreases 
with increasing Q 2 still a considerable fraction of all events is characterized by photons with 
rather large virtualities. These events cannot be expected to be described reliably within the 
present approach. In Figs. [7|a and 0c the energy dependence of the average total, positively 
and negatively charged hadron multiplicities are shown. For both target nuclei the data are 
well reproduced by the model. The multiplicities of charged hadrons are shown separately for 
the forward and backward region of the photon-nucleon cm. frame in Figs. 0b and [7|d. The 
calculated multiplicities in the backward region which are strongly affected by target associ- 
ated particle production and, therefore, by the laboratory momentum cut applied to the final 
state hadrons, are slightly higher than the data whereas those in forward direction are lower. 
Multiplicities averaged over all energies of the photon-nucleus interaction are compared to 
the corresponding experimental values in Table 0. 

3.2.3 Inclusive particle distributions 

In Fig. [| the model results for the pseudorapidity distributions of shower particles from 
muon-emulsion interactions at 150 GeV are compared to data again from measurements 
by Hand et al. ^7|. The kinematic range of the experiment is 0.6 < Q 2 < 21 GeV 2 and 



2.5 < W < 16.5 GeV. The distributions shown for muon-emulsion interaction cover different 
but overlapping VT-ranges: in a) data and MC results are plotted for 9 < W < 14 GeV and in 
b) for W > 10 GeV ((Q 2 ) = 2.7 GeV 2 ). Let us first compare the two data sets which include 



43 events in a) and 47 events in b) ||47|| . One obvious difference is the peak in the distribution 
in a) at 77 = 1.25 which is not present in b). The only energy range which is not covered by the 
data in Fig. ||a as compared to |]b is 9 < W < 10 GeV. Therefore, we assume that this peak is 



due to statistical uncertainties within the experiment (c.f. also our discussion in Sect . |3. 2. 2 ). 



Comparing the model results to the data we note that agreement within the statistical errors 
is obtained for r\ > 2. Disregarding the above mentioned peak, at lower pseudorapidities we 
underestimate the measured distributions slightly. Furthermore the comparison suggests that 
both, the measured and calculated distributions, agree in normalization (i.e. in the average 
multiplicity of shower particles) but the maximum of the calculated distribution appears at 
somewhat higher pseudorapidities than the measured one. In order to understand this dis- 
crepancy we compare in Fig. ||a in addition results for pion-emulsion interactions at 60 GeV 
to data ||4"9|. This energy corresponds to the average photon laboratory energy of the distri- 



bution for the muon projectile in this Figure. The cuts on both data sample are similar: only 
events with at least three heavy particles are considered and shower particles are defined by 
their Lorentz-/? value as (3 > 0.7. Our results for 7r~-emulsion interactions agree in shape 
and in position of the maximum to the data. Therefore and with respect to conclusions from 
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comparisons with data drawn further below we attribute the discrepancy in the position of 
the maximum to the statistical uncertainties of the experiment. 

Turning again to muon-deuterium and muon-xenon interactions at 490 GeV we compare 
in Figs. | and |T0| rapidity distributions of positively (a) and negatively charged hadrons (b) in 
the photon-nucleon cm. frame to data |4"8fl . The comparisons are shown for three ranges of the 
photon-nucleon cm. energy W. Taking into account that the treatment of high Q 2 -values in 
our model J15 might be too simplified, the description of the data by the model is satisfactory 
with the exception of the production of positively charged particles in the target fragmentation 
region of the muon-xenon interaction (Fig. [l0|a, y « —3). The peak in the calculated distri- 
butions clearly reflects the production of target associated particles by intranuclear cascade 
processes which are present in the distributions if a momentum cut as low as 200 MeV/c in the 
laboratory frame is applied to the results. The multiplicity seen in the target fragmentation 
region depends strongly on this cutoff. We assume that the differences might be due to addi- 
tional kinematic cuts applied to the data or due to experimental uncertainties (with respect 
to the momentum cutoff) for the following reasons: (i) the models for slow particle production 
implemented in Dtunuc 2.0 are in good agreement with data from hadron-nucleus as well as 
nucleus-nucleus collisions p2|, ^3|, (ii) the dependence of shadowing on the photon virtuality 
and energy is qualitatively understood within our model and describes corresponding data 
of the E665-Collab. reasonably well fli~0"| , and (iii) the rapidity distributions of negatively 
charged particles and of positively charged particles outside the target fragmentation region 
in muon-xenon interactions and of charged particles in muon-deuterium interactions agree 
with E665-data. 

Distributions of charged hadrons from muon-deuterium and muon-xenon interactions at 
490 GeV were measured as function of z = Ejv [E and v being the secondary hadron energy 
and the photon energy in the target rest frame, resp.) by Adams et al. |]50|| . In Figs. |TT| a and 



TT| b our results on the ^-distributions are shown together with these data. Both, model results 
and data, are restricted to the shadowing region, i.e. to XBj < 0.005 and Q 2 < 1 GeV 2 . We 
find a good agreement in the whole z-range. 

Feynman-x distributions are usually studied in terms of the one-particle inclusive variable 
F(xp) defined as 

F(x.) - 1 2E da x,- 2P|1 (10) 

F{xf) ~ XF ~ w (10) 

The quantities E, pn, and W = y/s denote the energy and longitudinal momentum in the 
photon-nucleon cm. system, and the photon-nucleon cm. energy, resp. The Feynman-x 
distribution of positively and negatively charged hadrons is given together with data [O] in Fig. 



TT| c for muon-deuterium interactions at 147 GeV. Here, the photon virtualities are restricted 
to the range 0.5 < Q 2 < 3 GeV 2 and the photon-nucleon cm. energies to W > 10 GeV. 

Finally, we compare in Fig. 0a the transverse momentum distributions of charged particles 
in different ranges of x? again to data on muon-deuterium interactions from Ref. ||51|| . Our 
results are presented as histograms whereas fits to measured ^-distributions are shown as 
continuous lines. The experimental uncertainties are increasing with p± and are at least of the 
order of the differences between the fits for positive and negative particles They do not 
allow to draw conclusions concerning the disagreement between model predictions and data 
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at large transverse momenta. The dependence of the average transverse momenta of charged 
particles on xp is shown in Fig. |l"2|b. 



3.2.4 The nuclear dependence of particle production 

In order to study the dependence of particle production on the mass number A of the nuclear 
target, inclusive single particle cross sections, like pseudorapidity (da/drj), transverse mo- 
mentum (da/dpY), or Feynman-x distributions (da/dxp), are usually fitted to a ^"-behavior. 
In the projectile fragmentation region of photon-nucleus collisions we expect a to approach 
unity (i) for small gg-nucleon cross sections, such as at low energies or for Q 2 > 0, (ii) at high 
energies for interactions becoming more point-like, and (iii) at large transverse momenta due 
to hard interactions. In the fragmentation region of the target nucleus a can exceed the value 
of one due to the formation zone intranuclear cascade, an effect being outside the scope of the 
present paper. Here, we restrict our discussion to the photon fragmentation region. 

In Fig. [IB] we show fits to results of the model on interactions of real and weakly virtual 
photons (Q 2 = 1 GeV 2 ) with carbon, aluminum, copper, silver, xenon, and gold. Only 
charged particles with Lorentz-/? > 0.7 (shower particles) are taken into consideration. The 
dependences of a on the pseudorapidity (a), on the transverse momentum (b), and on the 
Feynman-x variable (c) are plotted for a laboratory energy of 250 GeV. At all values of 
the studied variables a is larger for virtual than for real photon projectiles. This stronger 
A-dependence of the cross sections for virtual photons is due to the <5 2 -dependence of the 
effective gg-nucleon cross sections gvn (Eq.(§)). It demonstrates that the previously discussed 
property of the model, the Q 2 -dependence of shadowing (see Fig. |5]), is clearly visible also in 
inclusive particle distributions. It might be interesting to study this effect also experimentally. 



3.3 Jet photoproduction off nuclei at HERA-energies 

Studying jet production in direct and resolved photon events in jp interactions striking dif- 
ferences have been observed by the experiments |j[ |8|, ||. Therefore one can expect to observe 
similar effects when comparing point-like and resolved photon processes in photon-nucleus 
interactions. In the following, properties of hadronic jets produced in interactions of real 
photons with nuclei are studied at energies which might be reached with nuclear beams at 
HERA. Assuming nuclei with equal numbers of protons and neutrons the nuclear beam will 
have an energy per nucleon of 410 GeV ||46|| . With an average photon energy of about half of 
the electron energy (27.6 GeV) the photon-nucleon cm. energy will be about 150 GeV. 

Before discussing jet production let us first compare transverse momentum distributions 
of charged particles in proton- and photon-nucleus collisions. They very clearly reflect dif- 
ferences between interactions with proton and the photon projectiles which are responsible 
for characteristic features of jet production in photon-nucleus collisions. As shown in Fig. [14] 
the p^-distribution of charged particles for photon-projectiles exhibits a less rapid decrease 
and extends to higher transverse momenta than the one for hadron projectiles. This property 
was already discussed in [T3| comparing hadron-hadron, photon-hadron, and photon-photon 



interactions. The reasons are (i) the photon may interact in direct processes, and (ii) the 
photon PDFs are "harder" as compared to the proton PDFs. 
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As it has been shown in the analysis of data from HERA on photon-proton collisions 
the study of jet production may give evidence for the subdivision into direct and resolved 
photon interactions. Therefore it can be expected that differences between hadron-nucleus 
and photon-nucleus collisions and characteristic features of point-like photon interactions 
emerge most clearly in a study of jet production. In the following, all observables used refer 
to the photon-nucleon cm. system. Particle jets are defined according to the Snowmass- 
convention |53| . In the plane of pseudorapidity r\ and azimuthal angle a jet is defined as 



a collection of particles contained in a cone of radius R — J (Ar/) 2 + (A0) 2 = 1. The jet 

transverse energy E± is taken as the sum of the particle transverse energies inside the cone. 
The jet pseudorapidity rj }et is calculated as ^-weighted average over the pseudorapidities of 
the particles belonging to the jet. 

In Fig. |T5| a we compare the transverse energy distributions of jets from proton-carbon 
and photon-carbon interactions. Similar to proton/photon-proton interactions [T3j, the 



distribution extends to higher transverse energies for photon-nucleus collisions as compared 
to proton-nucleus collisions. Differences between the two channels are even more pronounced 
within the pseudorapidity distributions of the jet axes. This is shown in Fig. |T5| b for jets 
with transverse energies higher than 6 GeV and 10 GeV. As initially mentioned, since the 
PDFs of the photon are harder than the PDFs of the proton and due to the point-like photon 
interactions we obtain in photon-nucleus collisions considerably more jets in forward direction 
(rf et > 0). In backward direction the soft underlying event, being more pronounced for proton 
projectiles, is responsible for a higher jet rate. 

The kinematic regions of jets from direct and resolved photon interactions might be sepa- 
rated if jet production is analysed as function of 



2E~ 



x obs = : — ^- — . (ll) 



E 7 denotes the photon energy, x can be considered as an estimate for the fraction of the pro- 
jectile momentum entering the hard scattering [§, In Fig. [I6| a we show x obs -distributions 
for proton-carbon, photon-carbon, and photon-sulfur interactions. Point-like photon interac- 
tions contribute exclusively to the region x obs > 0.6. The specific properties of jets and of the 
underlying event are usually studied in terms of jet profiles. This is shown for photon-carbon 
interactions and different x obs -bins in Fig. |T6"| b where we plot the average transverse energy as 
a function of the distance in pseudorapidity from the jet axis. Since in direct photon inter- 
actions there is no photon remnant which could contribute to the jet pedestal, the transverse 
energy outside the jets is decreasing with increasing x obs . Of course, as the mass number 
of the target nucleus increases the soft underlying event becomes more pronounced, i.e. the 
average transverse energy of the jet pedestal increases. This is demonstrated in Fig. |T6| c for 
four different target nuclei. Here, all jets with E± > 6 GeV are included. 



4 Summary and conclusions 

The two-component Dual Parton Model is applied to multiparticle photoproduction off nu- 
clei. By incorporating the PHOJET event generator for hadron-nucleon and photon-nucleon 
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interactions, it is now possible to describe with the event generator Dtunuc 2.0 consistently 
cross sections and particle production in high energy hadron-nucleus, real and weakly virtual 
photon-nucleus, and nucleus-nucleus collisions (the latter are not discussed in this paper). 
In |g it has been shown that the model correctly describes the shadowing behavior of high 



energy photon-nucleus cross sections. Based on this fact and on the overall good agreement 
of model results and data in all main aspects of particle production in photon-hadron and 
hadron-nucleus collisions no further freedom exists in the model for its application to multi- 
particle photoproduction off nuclei. It is therefore a severe test of the model. Unfortunately, 
only a few data are available in the kinematic region Q 2 < s, Q 2 < 9 GeV 2 to which the 
predictions of Dtunuc 2.0 can be compared. As discussed in this paper, model results are 
qualitatively consistent with many features of these data. 

Within the model the photon is treated (i) as a resolved gg-state interacting with target 
nucleons according to the GVDM and the Gribov-Glauber approximation in multiple soft 
and hard scattering processes and (ii) as a point-like object interacting in a single hard scat- 
tering with one target nucleon (anomalous component of the photon PDF and direct photon 
interactions). On the basis of this treatment we obtain the following results: 

• As expected from studies of photon-nucleus cross sections, also particle production off 
nuclei clearly shows decreasing shadowing with increasing photon virtualities. 

• At energies of present fixed target experiments, inclusive single particle cross sections 
become proportional to A a ,a ~ 1, already at Q 2 > 1 GeV 2 due to decreasing qq- 
nucleon cross sections with increasing Q 2 . However, more data on particle production in 
photon-nucleus collisions would be needed for a detailed investigation of the transition 
region. 

• In analogy to observations in photon-proton collisions at HERA we argue that also 
in photoproduction off nuclei the dual nature of the photon and, therefore, point-like 
photon interactions show up most clearly in hadronic jet production. As examples 
we present transverse energy and pseudorapidity distributions of jets and we study jet 
profiles in real photon-nucleus collisions at energies which might be reached with nuclear 
beams at HERA. The strong dependence of the jet profiles and of the underlying event 
on the fraction of the momentum of the photon going into jets which was first observed 
at HERA, can be expected to be present also in photon-nucleus collisions. 
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Tables 



Table 1: Multiplicities of shower particles {(3 >0.7) in interactions of pions and protons with 
carbon-, copper-, and lead-nuclei at 50, 100, and 150 GeV are compared to measurements |43| . 





£ Lab =50 GeV 


£ La b=100 GeV 


S Lab =150 GeV 




Dtunuc 2.0 


Exp. 


Dtunuc 2.0 


Exp. 


Dtunuc 2.0 


Exp. 


7T+-C 


7.3 


7.62±0.14 


9.0 


9.19±0.17 


9.9 


10.01±0.18 


7T+-CU 


8.6 


8.81±0.23 


10.5 


10.41±0.27 


11.8 


11.57±0.30 


7T+-Pb 


10.1 


10.11±0.38 


12.2 


11.54±0.44 


14.3 


13.08±0.50 


p-C 


7.7 


7.88±0.15 


9.2 


9.25±0.18 


10.3 


10.58±0.19 


p-Cu 


9.8 


9.52±0.24 


11.9 


11.40±0.29 


13.3 


12.93±0.33 


p-Pb 


11.1 


11.31±0.43 


13.8 




15.9 


14.95±0.57 
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Table 2: Results on average multiplicities of all charged (-/V c h), positive (N + ), negative (iV_), 
forward (Np), and backward {Nb) hadrons (p > 200 MeV/c) for muon-deuterium and muon- 
xenon interactions at 490 GeV are compared to data of the E665-Collab. [48 . 







D 2 




Xe 




Dtunuc 2.0 


Exp. 


Dtunuc 2.0 


Exp. 


(N ch ) 


7.7 


7.83±0.07 


10.8 


9.99±0.13 


(N + ) 


4.1 


4.16±0.05 


6.8 


6.00±0.11 


(N-) 


3.6 


3.67±0.05 


4.0 


3.99±0.06 


(N F ) 


4.1 


4.60±0.05 


4.2 


4.70±0.07 


(Nb) 


3.6 


3.23±0.04 


6.6 


5.29±0.11 
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Figure Captions 



1. Color flow picture of a single reggeon exchange graph (a) and the corresponding unitarity 
cut (b). In c) and d) the same is shown for a single pomeron exchange. 



2. Comparison of model results to data on pion-emulsion interactions [44]. In a) shower 
particle multiplicity distributions are shown for a laboratory energy of 525 GeV. In 
b) pseudorapidity distributions of shower particles are given for two different energies: 
60 GeV and 525 GeV. 

3. Model predictions on transverse momentum distributions of 7r°'s in proton-gold inter- 
actions at 200 GeV are compared to data of the WA80-Collab. f45 |. 



4. Average numbers of copper nucleons (ut) struck by projectile pions, real and virtual 
(Q 2 =2 GeV 2 ) photons are given as function of the projectile laboratory energy (a). In 
b) and c) average multiplicities of shower (N s ) and heavy particles (N h ) are presented. 
Shower particles are defined as singly charged particles with j3 > 0.7. Heavy particles 
are charged particles with (3 < 0.7 except residual copper nuclei. 

5. In a) the Q 2 -dependence of the average number of copper nucleons interacting with the 
photon (i/ t ) is given for a photon-nucleon cm. energy of 150 GeV. In b) the dependence 
of the average heavy particle multiplicity on v t is shown for proton, pion, real photon, 
and weakly virtual photon (Q 2 = 2 GeV 2 ) projectiles at the same energy. 

6. Average numbers of charged hadrons (/3 > 0.7) in muon-emulsion interaction events at 
150 GeV with more than two heavy particles in the final state are given as function of 



the inverse of the Bjorken-x variable and are compared to data [47 . 



7. The dependence of average multiplicities of charged hadrons on the squared energy of the 
virtual photon-nucleon cm. system is shown for muon-deuterium (a,b) and muon-xenon 
interactions (c,d) at 490 GeV and is compared to measurements of the E665-Collab. ||48|| . 
In a) and c) we give the average multiplicities of all (N ch ), of positively (N + ), and of 
negatively (AL) charged hadrons. In b) and d) the average multiplicities of all charged 
hadrons are shown for positive (Np) and negative (Nb) cm. rapidities. 

8. Pseudorapidity distributions of charged hadrons {(3 > 0.7) from muon-emulsion inter- 
actions at 150 GeV are shown together with data |[47|| . In a) the photon-nucleon cm. 
energy range is restricted to 9 < W < 14 GeV. In addition, the pseudorapidity dis- 
tribution of charged hadrons from pion-emulsion interactions at 60 GeV is plotted and 



compared to data | 49|j . In b) the results for W > 10 GeV are given. In both distributions 



only events with Nh > 3 are included. 

9. Rapidity distributions of positive (a) and negative hadrons (b) from muon-deuterium 
interactions at 490 GeV are shown together with data HSfl for different ranges of the 
virtual photon-nucleon cm. energy W. 

10. As in Fig. here for muon-xenon interactions. 
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11. Energy- (z-) distributions of charged hadrons from muon-deuterium (a) and muon- 



xenon interactions (b) at 490 GeV are compared to measurements of the E665-Collab. | 50fl . 
Corresponding to the experimental cuts applied to the data the MC-results are restricted 
to XBj < 0.005, Q 2 < 1 GeV 2 , and v > 100 GeV. In c) model results for Feynman-x 
distributions of charged hadrons from muon-deuterium interactions at -Exab = 147 GeV 



are compared to data [51 



12. In a) transverse momentum distributions of charged hadrons from muon-deuterium 
interactions at 147 GeV are presented for different Feynman-x ranges. The model results 



(histograms) are plotted together with fits to data as given in [51]. In b) the x-p- 
dependence of the average transverse momenta of charged hadrons are compared to 
data Ml. 



13. The power a of the nuclear dependence of inclusive charged particle cross sections is 
shown as function of the pseudorapidity (a), the transverse momentum (b), and the 
Feynman-x variable (c) for real photon and weakly virtual photon projectiles of 250 GeV 
laboratory energy. 

14. Transverse momentum distributions of charged hadrons from proton- and photon- 
carbon interactions at a proton/photon-nucleon cm. energy of 150 GeV. 

15. In a) transverse energy distributions are presented for hadronic jets from proton-carbon 
and real photon-carbon interactions at an energy of 150 GeV in the proton/photon- 
nucleon cm. system. In b) the pseudorapidity distributions of the jet axes are shown 
for two lower cuts in transverse energy. 

16. Dependence of jet production on the fraction of the projectile momentum carried by 
the jets (x obs ) at a proton/photon-nucleon cm. energy of 150 GeV: In a) the x obs - 
distributions of proton- and real photon-carbon and of real photon-sulfur interactions 
are plotted. Jet profiles as function of the distance in pseudorapidity from the jet axis 
are presented in b) for different x obs -bins. In c) the jet profiles are shown for different 
target nuclei and all x obs - values. 
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